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An  electron  accelerati’n  method  i'l  investigated  which  employs  a  short  (t,.  ~  2-k  '  ~  I  picosec).  high 

p<iwer  |P>  lo''’  W),  single  frequency  laser  pulse  to  generate  large  amplitude  (E  >  1  GeV/ml  plasma  waves 
(Wakefields).  At  sufficiently  high  laser  powers  (P  >  17  fu.-/,,-,  *'  OW),  relativistic  optical  guiding  may  be 
used  to  prevent  the  pulse  from  diffracting  within  the  plasma.  / 
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LASER  WAKEFIELD  ACCELERATION  AND 
RELATIVISTIC  OPTICAL  GUIDING 


Introduction 

It  has  been  suggested  that  the  next  generation  of  high  energy  linear 
electron  accelerators  utilize  the  extremely  high  gradients  associated  with 
plasma  waves.  Excitation  of  plasma  waves  having  gradients  as  high  as 
several  tens  of  flpV/m  can  be  accomplished  in  a  number  cf  ..elated  ways.  The 
plasma  based  acceleration  schemes  which  have  received  the  most  attention 
are  the  plasma  beat  wave  accelerator^  (PBWA)  and  the  plasma  wakefield 

9 

accelerator  (PVFA). 

The  purpose  of  this  report  is  to  propose  a  laser  plasma  electron 
3 

acceleration  scheme  which  utilizes  a  relativistic  optical  guiding 

mechanism.  Relativistic  optical  guiding"^’^  may  allow  a  sufficiently  high 

power  laser  pulse  to  propagate  long  distances  within  a  plasma.  The 

principle  of  this  optically  guided  laser  wakefield  accelerator  (LVFA)  is 

that  a  short  (x,  ~  2ii/(a)  ~  1  picosec),  high  power  (P  >  10^^  W),  single 

up 

frequency  laser  pulse  could  propagate  long  distances  in  a  plasma  and 
produce  accelerating  Wakefields  in  a  manner  analogous  to  that  in  the  PWFA 
(see  Fig.  1).  In  the  LVFA,  however,  the  plasma  responds  to  tiie 
ponderomoti ve  forces  of  the  laser  pulse  as  opposed  to  the  self-fields  of 
the  electron  beam  as  in  the  PWFA.  In  addition,  in  the  LVFA  the  plasma  wave 
is  not  resonantly  excited  as  it  is  in  the  PBWA.  Therefore,  the  plasma 
density  in  the  LWFA  concept  does  not  have  to  be  finely  adjusted  to  achieve 
large  amplitude  accelerating  fields.  The  idea  of  generating  a  plasma  wave 
using  a  single  frequency,  short  pulse  laser  was  suggested  by  Tajima  and 
Dawson,^  but  apparently  was  not  pursued.  Mote  detailed  cons i dci a t i on  of 
the  laser  propagation  issues,  along  with  recent  advances  in  laser 
technology,  indicate  that  the  single  frequency,  slioi  t  pulse  L.WFA  logethei 
with  relativistic  optical  guiding  may  have  advani.ages  over  the  PBWA  and 
PWFA  schemes. 

Manuscript  approved  July  1988 
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In  the  PBUA  the  plasma  wave  is  excited  by  the  beating  of  two 
relatively  low  power,  long  pulse  laser  beams  having  a  frequency  difference 
equal  to  the  plasma  frequency.  The  beat  (ponderomotive)  wave  resonantly 
drives  the  plasma  wave  Lo  large  amplitudes.  In  the  PWFA  concept,  a  low 
energy,  high  current  relativistic  electron  beam  (driver)  having  an 
appropriate  current  profile  travels  through  a  plasma  leaving  behind  a  large 
amplitude  plasma  wave  (wakefield).  The  wakefield  accelerates  a  second  low 
current,  high  energy  relativistic  electron  beam.  A  necessary  criteria  for 
successful  operation  of  either  tlie  PBWA  or  the  PWFA  is  that  the  driver, 
i.e.,  radiation  or  electron  beam,  must  be  capable  of  propagating  a 
sufficiently  long  distance  within  the  plasma. 

Both  the  PBWA  and  the  PWFA  concepts  have  a  number  of  unresolved 

issues.  In  the  PBWA,  these  include  fine  tuning  of  the  laser  frequencies 

and  plasma  density  to  within  a  fraction  of  a  percent  to  allow  for  resonant 

growth  of  the  plasma  wavc.^  Also  the  laser  beams  must  propagate  large 

distances  within  the  plasma,  avoiding  i)  diffraction,  ii)  laser-plasma 

instabilities,  iii)  phase  detuning  between  the  plasma  waves  and  the 

accelerated  electrons,  as  well  as  iv)  energy  depletion  of  the  driving  laser 
7 

beams.  The  problems  regarding  the  PWFA  involve  the  technology  of 
producing  a  high  current  driving  beam  with  a  slow  rise  time  and  a  very 

g 

rapid  fall  time,  of  the  order  of  picoseconds,  as  well  as  the  stable 
propagation  of  such  a  beam  over  large  distances  within  the  plasma. 

Multiple  acceleration  stages,  all  sequentially  phase  synchronized,  have 
been  proposed  to  overcome  the  propagation  distance  limitation  in  both  the 
PBWA  and  PWFA.  Mul t i -s tagi ng  appears  to  be  extremely  difficult  from  a 
practical  point  of  view. 
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optical  Guiding 

The  need  for  optical  guiding  in  the  LVFA  becomes  apparent  when  the 
various  limitations  placed  on  the  acceleration  distance  are  considered. 
One  limitation  on  the  acceleration  distance  is  the  diffraction  length, 
L^,  which  characterizes  the  distance  over  which  the  laser  beam  spreads 
transversely.  In  the  absence  of  some  form  of  optical  guiding,  tiic 

2 

diffraction  length  is  given  by  the  vacuum  Rayleigh  length,  =  nr^^/X, 


1,7 


where  rj^  is  the  laser  spot  size  and  X  is  the  wavelength.  Another 

limitation  or  the  acceleration  length  is  the  phase  detuning  distance, 

L  =  2y?X  =  2(l  +  X^/4r^)  \x  /X)^X  ,  where  y.^  -  (l-v^/c^),  v  is  the  group 
t  Lp  'pL  'p'p  L  g  g 

velocity  of  the  laser  pulse  and  X^  is  the  plasma  wavelength.  The  phase 
detuning  length  is  the  distance  over  which  an  ultra  relativistic  electron 
outruns  the  Wakefield  of  the  radiation  pulse  and  no  longer  gains  energy. 

9 

In  addition  to  and  there  is  also  the  laser  depletion  length, 

I,  =  =  C,(C,/X)^/a^  ,  where  E,  is  the  laser  electric  field,  [,  is 

pLLzL  L'Lo  L  ’L 

the  laser  pulse  length,  is  the  axial  wake  electric  field  and  aj^^  is  the 
normalized  vector  potential  amplitude  of  the  radiation  field. 


|e |A.  /(m  c  ).  When  the  pulse  travels  a  distance  L  ,  the  energy 


in 


'Lo  Lo  '  o  '  p’ 

the  trailing  plasma  wakefield  becomes  comparable  to  the  laser  pulse  energy. 

Typical  values  for  ,  L^.  and  are  ~  1  m,  ~  100  m  and  ~  1000  m 

respectively.  In  obtaining  these  estimates  the  following  parameters  wete 


used : 


1  pm ,  a 


Lo 


0.5  and  C, 


X^  -0.5  mm.  The  primary 


limitation  on  the  acceleration  distance  is  due  to  diffraction.  L,. 

d 

Clearly,  some  form  of  optical  guiding  within  tlie  plasma  is  necessary  to 
avoid  the  need  for  multi  stage  acceleration. 

The  optical  guiding  mechanism  which  may  be  appiopi late  for  the 
intense,  short  laser  pulse  in  tire  LUFA  is  that  of  relativistic  guiding.^'  ^ 
Physically,  relativistic  guiding  lesults  from  the  quiver  motion  of  the 


3 


plasma  electrons  in  the  radiation  field,  v  =  ca, /y,,  where  Y,(r)  = 

Q  ij  X  X 

2  1/2 

(l+aj^(r))  .  This  gives  an  index  of  refraction  n(r)  = 

2  2  1/2 

(l-fWp^/w  )/Yj^(r))  ,  where  is  the  ambient  electron  plasma  frequency 

and  w  is  the  laser  frequency.  If  the  radiation  beam  is  peaked  on  axis, 

then  3n/8r  <  0,  which  is  a  necessary  requirement  for  refractive  guiding  to 

occur.  Relativistic  optical  guiding  occurs  on  a  fast  time  scale  of  order 

w  S  hence,  it  can  affect  short  pulse  radiation,  <<  e  /c  < 

L  ~  p 

Using  the  ray  equations  from  geometric  optics,  it  is  possible  to 

derive  an  envelope  equation^  for  the  evolution  of  the  normalized  spot  size 

X  s  of  tho  radiation  beam,  where  r^^^  is  the  initial  spot  size. 

The  envelope  equation  is  of  the  form  of  a  particle  moving  in  an  effective 
2  2 

potential,  d  x/dt  =  -V^  3V/3x.  The  effective  potential  V(x)  is  given 

by  3V/3x  +  16o(x(g(x)  -  2  ln(g(x)/2  +  l )  J ,  where 

=  (2c2/(a)r2^aJ^))2,  «  =  ,  and  g(x)  =  (Ux-^)^^^-!. 

Analysis^  indicates  that  the  effective  potential  contains  a  minimum 

provided  «  >  1,  thus  allowing  for  matched  beam  (constant  spot  size) 

solutions.  Physically,  a  can  be  written,  in  terms  of  the  laser  power  P,  as 

2 

a  =  P/P  ,  where  P  =  17(w/w  )  GW  is  the  critical  power  threshold  for 
c  r  c  r  p 

relativistic  optical  guiding.  The  high  power  levels  needed  for 

relativistic  optical  guiding  in  plasmas  are  consistent  with  the  intense 

laser  pulses  needed  in  the  LWFA. 

Two  points  should  be  mentioned  with  regard  to  the  propagation  of 
finite  length  pulses  of  duration  f^/c  <  •  T'’’®  first  is  that 

relativistic  optical  guiding  may  also  lead  to  "pulse  clipping".  That  is, 

the  front  and  back  regions  of  the  pulse  where  P  <  P^ ^  will  not  be  guided 


but  instead  will  diffract  away,  leaving  a  slini 


l  I  t  VJ  1 


St.  Only  the 


central  region  of  the  pulse,  where  P  >  j.i  will  propagate.  The  second 
point  concerns  longitudinal  dispersive  spreading.  It  can  be  shown  that 
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after  propagating  a  detuning  length  L^,  the  intrinsic  frequency  spread  of 
the  beam  Aw  causes  the  pulse  to  spread  by  the  amount  AC,  =  2(Aw/w)X  . 

L  p 

Since  |Aw/w|  <<  1,  longitudinal  dispersive  spreading  should  not  be  a 
problem. 

Acceleration  Mechanism 

In  the  relativistically  guided  LWFA  concept  the  short  pulse,  high 
power  laser  beam  provides  both  a  radial  and  axial  ponderomot i ve  force  on 
the  plasma  electrons.  The  radial  ponderomotive  force  expels  electrons 
radially  outward  while  the  front  (back)  of  the  laser  pulse  exerts  a  forward 
(backward)  force  on  the  electrons.  In  this  sense,  the  laser  pulse  acts 
approximately  like  a  negatively  charged  macro  particle  propagating  through 
the  plasma  (see  Fig.  1).  As  the  plasma  electrons  flow  around  the  laser 
pulse,  large  amplitude  plasma  waves  are  generated. 

The  ponderomotive  force,  exerted  by  the  laser  pulse  on  the  plasma, 

moves  at  the  pulse's  group  velocity  and  is  given  by  =  je  |  V'J'j^(  r ,  7 .  t ) , 

2  2 

wlieie  the  ponderomotive  potential  is  =  -m^c  aj^/(2  |e  | ) .  Note  that  tlie 
axial  ponderomotive  force  from  the  laser  pulse  cannot  be  used  directly  to 
accelerate  electrons  to  high  energies.  The  ponderomotive  force  on  the 
accelerated  electrons  is  smaller  than  that  on  the  plasma  electrons  by  the 
factor  1/y»  where  y  is  the  relativistic  factor  associated  with  the 
accelerated  electrons.  The  laser  pulse  must  first  excite  a  plasma  wave 
which,  in  turn,  can  be  used  for  acceleration.  In  this  analysis  the  laser 
beam  is  assumed  to  be  circularly  polarixed,  although  a  linearly  polarized 
laser,  apart  from  generating  harmonics,  would  liave  been  equally 
satisfactory. 
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The  wave  equation  for  the  plasma  response  or  wakefield  is 


,  ,  3^E 

'E  -  Tri- “> 

c  3t  c  ^ 


where  SJp  and  Sn^  are  the  plasma  response  electron  current  and  number 

density  respectively.  It  proves  convenient  to  perform  an  algebraic 

transformation  to  the  speed  of  light  frame  (C  =  z-ct,  x  =  t).  The 

transformation  should  actually  be  to  the  laser  pulse  group  velocity  frame, 

but  the  differences  can  be  neglected  for  the  present  purposes. 

Furthermore,  a  temporal  steady  state,  9/9t  =  0,  in  the  laser  pulse  frame 

2 

is  assumed.  It  can  be  shown  that  for  short  laser  pulses  with  a,  /2  <<  1, 

Lo 

the  plasma  quantities  remain  linear  and  nonrela t i vis t ic .  The  plasma, 
therefore,  is  assumed  to  be  described  by  the  linear,  nonrelativistic,  cold 
fluid  equations.  Using  this  fluid  response,  SJ^  and  Sn^,  in  the  wave 
equation,  the  plasma  response  fields  and  density  are  given  by 


E(r,C) 


5n^(r,C)  =  - 


e  n 


PO 


2 

m  c 
o 


V^<J>j(r,C), 


(2) 


(3) 


where  k  =  w  /c  and  V  =  e  9/3  +  e  9/9C.  Note  that  even  in  the  two- 

p  po  ~  r  r  z 

dimensional  case  the  response  field,  E,  is  derivable  from  a  scalar 
potential  and  hence,  there  is  no  response  magnetic  field. 

From  (2),  the  axial  wakefield  is  given  by 


E^(r,0 

I 


cr» 


J 


kp(C  C)<t>,  (r,C'  )dC'  . 


(4) 


From  (2)  and  (3)  it  can  be  shown  that  the  transverse  wakefield  and 

- 1  2 

density  are  given  by  9E^/9C  "  9E^/9r  and  9Snp/9C  -  (4n.|e|)  V  E^. 


plasma 
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As  an  illustration,  consider  a  laser  pulse  profile  of  the  form 

2  2 

aj^(r,C)  =  aj^^sin(iiC/?j^)exp(-r  0  £  C  <  and  0  otherwise.  Then 

the  axial  Wakefield  and  response  plasma  density  within  the  laser  pulse, 

0  <  C  <  and  behind  the  pulse,  C  <  0,  are  given  by 


2ii  k  (|>,  (r) 

E  (C,r)  =  -  - 

4n  -k  ef 
P  L 


^sin  kpCj^(l-C)  +  h  sin  ^k^C j^C/hj  ,  (5) 


6n  (z)  2  e  it'^i,  (r) 

p  Lo 


m  c  4rt  -k  C,  1 
o  t  p  L  ; 


P  L 


cos  kpCj^(l-0 


-  1  -  h2(cos(k  C^C/h 


-  1  -(6) 


where  Z  =  C/C^,  k^  =  2n/Xp,  ^  ^^Lo  '  2r^/r^)  and 

where  h  =  k  C, /2n  for  0  <  C  <  1  and  h  =  1  for  C  <  0.  The  transverse 

Wakefield  is  easily  calculated  from  (5)  by  the  relation  3E^/RC  =  8E^/0r. 

It  can  be  shown,  as  is  true  of  PUFA,  that  there  exists  a  region  of  length 

X  /4  in  the  laser  pulse  frame  over  which  the  accelerated  elections 
P 

experience  both  an  accelerating  axial  field  as  well  as  a  focusing  radial 


field. 


The  axial  Wakefield  in  (5)  is  maximum  when  the  laser  pulse  length  is 


nearly  equal  to  the  plasma  wavelengtii,  '  ^  =  X^.  lor  <  ^  X^,  the  rnaximiii'; 

acceleiating  field  is  approximately  n  times  larger  than  the  maximum 

9 

ponderomot  i  ve  axial  field  E  =  nF.  ,  -  tt‘ 41,  .  It  can  be  slio  1; 

^  z,max  pond, max  Lo  L 

that  the  maximum  accelerating  field  is  fairly  insensitive  to  changes  in  the 
laser  pulse  length  and/or  the  ambient  plasma  density.  It  stiould  he  noted 
that  (5)  and  (6)  also  indicate  that  it  is  possible  to  opeiate  the  lUFA  in  a, 
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wakeless"  regime  (i.e.,  the  plasma  response  is  nonzero  only  within  the 


region  of  the  laser  pulse)  when  where  m  is  an  integer  >  2. 


Numerical  Results 

The  results  for  the  plasma  response  given  by  (5)  and  (6)  are  plotted 

2 

in  Fig.  2  for  the  parameters  ^  "  0.03  cm,  a^^  =  0.31  and 

=  0.038  cm.  The  values  of  a^^  and  r^^  are  those  required^  for  a 
relativistic  optical  guided  beam  when  a  =  =  ^•2-  The  axial  wakefield 

is  shown  by  the  solid  curve  and  the  density  wake  is  shown  by  the  dashed 
curve.  The  maximum  accelerating  gradient  for  this  example  is  2.6  GeV/m. 
Recall  that  the  laser  pulse  extends  over  the  region  0  <  C  <  1. 

In  order  to  further  examine  the  principles  of  the  LWFA,  a  full  scale 
simulation  was  performed  using  the  electromagnetic  particle  code^*^  FRIEZR. 
FRIEZR  is  a  2  1/2D,  fully  relativistic,  electromagnetic  PIC  code  for 
electrons  with  a  fluid  ion  background.  The  simulation  is  carried  out  in 
tlie  transtormed  laboratory  frame  of  C  -  z-ct.  The  laser  field  was  modeled 
by  a  fixed  external  ponderomot i ve  force  moving  at  the  speed  of  light.  The 
resulting  axial  wakefield  is  shown  in  Fig.  1  for  the  same  parameters  as 
used  in  Fig.  2.  The  results  shown  in  Fig.  3  are  in  good  agreement  witli 
analytic  theoiy. 


D i scuss i on 

The  above  analysis  indicates  that  the  I.VFA  is  capable  of  generating 
acceletation  giadients  on  the  order  of  a  few  GeV/m  !)y  [)i  opaga  t  i  n.g  a  single, 
siioi  t  pulse,  liigh  pn<,'ei  lasei  beam  ihiough  a  plastr'a.  Equation  (b)  gives  a 

0  I 

maximum  acceletation  giadient  of  E  -m  (etta,  )'  /(2  e  i  ,  ).  In  addititjii. 
relativistic  optical  guiding  occuis  fot  ,su  f  f  i  c  i  eti  f  1  y  high  ladiafion  [toweis, 
P  >  ^ .  If  the  ladiation  pulse  is  rqitically  guided,  the  acceletation 
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distance  will  be  limited  to  the  phase  detuning  length,  L^,  instead  of  the 

much  shorter  free  space  Rayleigh  length,  L^.  This  indicates  a  maximum 

single  stage  energy  gain  of  S=  ~  ^“^p  “  =  ^'^^cr' 

Table  1  summarizes  these  results  for  a  CO2,  an  Nd  glass  and  a  KrF  laser, 

each  of  1  psec  pulse  duration.  In  each  case  a  =  1.2  which  implies 
2 

=  0.31  and  r^  =  0.038  cm  for  a  matched  beam  propagation  in  the 

relativistic  optically  guided^  propagation  mode. 

The  present  analysis  of  relativistic  optical  guiding  neglects  the 

effects  of  the  electron  density  response  on  the  laser  pulse.  Such  an 

2 

approximation  is  appropriate  when  Sn^/n^^  <<  a^^/2.  For  the  present 

analysis,  however,  this  condition  is  only  marginally  satisfied  for 

parameters  of  interest.  In  addition,  laser  plasma  instabilities  '  such 

as  the  f i lamen ta t ion ,  self-modulat ion  or  Raman  scatteiing  processes  have 

not  bi.en  considered  for  relat  ivist  ically  guided  sliort  pulses.  It  is 

anticipated  that  by  keeping  the  dimensions  of  the  laser  pulse  small, 

~  r^  <  X^,  the  effects  of  tliese  instabilities  may  be  minimized.  For 

1 2 

example,  Raman  scattering  processes  occur  througli  tlie  development  of 

plasma  waves  within  the  laser  pulse.  Since  the  length  scale  for  tlie 

development  of  plasma  waves  is  X^,  such  effects  may  lie  suppressed  if 

C,  <  X  .  In  addition,  relativistic  f i lamen ta t i on ^ ^  is  a  result  of  unstaliJe 

transverse  modes  with  >  X^.  Again,  this  instability  may  be  suppressed 

in  laser  pulses  with  r,  ~  X  .  Furtheimore,  landom  fluctuations  in  the 

L  p 

plasma  density  will  result  in  spreading  of  the  laser  spot  size.  A  mote 
se 1 f  -  cons i s t en t  model  of  relativistic  optical  guiding  tor  finite  pulse 
lengths  is  curiently  being  pursued  by  tlu?  authors. 

The  I.WFA  may  have  advantages  ovei  both  the  PWFA  and  tlu'  PBWA.  Fiu 
example,  in  the  PVFA,  it  is  necessary  to  use  a  high  cuiieni  (tens  of  kA) 
driving  electron  beam  with  a  long  lise  time  (>>  and  a  rapid  fall  time 
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plagued  by  the  usual  laser-plasma  instabilities.  In  addition,  such 
resonant  amplification  requires  fine  tuning  between  the  frequency 
differences  of  the  two  lasers  and  the  plasma  frequency.^  This  fine  tuning 
which  is  not  necessary  in  the  LUFA,  may  be  difficult  to  achieve  in 


practice.  Although  the  maximum  gradients  attainable  in  the  LWFA  may  be 
lower  than  in  the  PBWA,  the  many  apparent  advantages  (i.e.,  relativistic 
optical  guiding,  stability  and  simplicity)  of  using  a  single,  intense, 
short  pulse  laser  beam,  makes  the  LUFA  an  attractive  acceleration  scheme. 
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Table  1  The  laser  power  (P),  diffraction  length  (free  space  Rayleigh 


length,  Ij),  detuning  length  (L^)  and  single  stage  energy  gain 

i^=  for  three  lasers:  CO2,  Nd  glass  and  KrF.  The 

parameters  are  chosen  to  correspond  to  a  relativistic  optically 

guided  beam  with  P/P  =  1.2,  a^  ^  0.31  and  r, -  0.038  cm.  Thi 
®  cr  Lo  L 

gives  an  acceleration  gradient  of  E  =2.6  GeV/m  for  C,  =  X  = 

®  ^  z  L  p 

0.03  cm. 
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both  the  radial  and  axial  directions. 


(GeV/m 


Fig.  3  The  axial  wakefield  obtained  from  the  particle  code  FRIEZR  for  the 
same  parameters  as  in  Fig.  2.  The  laser  pulse  is  modeled  by  a 
fixed  ponderomot ive  force  which  extends  over  the  region 
0  <  (z-ct)/e^  <  1. 
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